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ABSTRACT 


The  present  paper  presents  an  analytical  study  of  the 
steady- state  velocity  field  induced  in  a  liquid  by  a  rotating 
cone.  As  a  consequence  of  the  experimental  observations, 
velocity  functions  which  very  accurately  define  the  flow  have 
been  developed.  The  functions  are  of  such  a  nature  that  the 
centerline  of  a  vortex  and  the  flow  about  the  centerline  are 
accurately  described.  The  boundary  conditions  at  the  surface 
of  the  rotating  cone  and  at  the  fixed  boundaries  are  satisfied. 
Also  the  continuity  equation  is  satisfied. 

Finally,  the  stress  field  in  the  fluid  has  been  studied 
in  terms  of  the  velocity  functions.  In  particular,  the  rela¬ 
tion  between  angular  velocity  of  cone  and  the  total  torque  on 
the  cone  was  determined  by  precise  experiments.  The  relation 
is  linear  and  a  discussion  of  the  implications  of  this  fact 
in  terms  of  the  nature  of  the  constitutive  equations  for  the 
liquid  is  presented.  The  force  of  the  evidence  is  that  the 
liquids  studied  are  Newtonian  and  the  viscosity  coefficients 
determined  by  experiment  agree  with  those  cited  by  various 


authorities. 
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I.  INTRODUCTION 

In  1953,  there  was  introduced  an  experimental  apparatus 
which  has  become  known  as  the  Ferranti-cone-plate  viscometer'*'. 
Determination  of  coefficients  of  viscosity  of  liquids  tested 
with  this  apparatus  is  based  upon  the  assumption  that  a  liquid 
always  flows  in  concentric  circular  paths  about  the  axis  of 
rotation  of  the  cone.  In  October  1962,  in  a  technical  paper 
entitled,  "A  Rotational  Fluid  Flow  Generator  for  Studies  In 

p 

Rheology"  ,  it  was  shown  that  actually  a  more  complex  flow 
pattern  exists  for  all  speeds  investigated.  It  was  further 
shown  that  the  observed  flow  is  stable  at  all  velocities  in 
the  range  studied,  even  down  to  zero  velocity.  Hence,  in 
particular  it  is  not  of  the  Taylor  instability  type  which  is 
generated  between  two  rotating  coaxial  cylinders. 

Motivated  by  these  discoveries,  It  was  decided  to  study 
the  actual  flow  pattern  and  to  make  a  comprehensive  analysis 
of  the  corresponding  steady-state  velocity  field.  An  effective 
method  available  for  the  study  is  based  upon  direct  observation 


of  the  flow  in  the  above  mentioned  generator.  Flow  visualiza¬ 
tion  techniques  and  torque  measurements  reveal  a  considerable 
amount  of  information  concerning  the  nature  of  the  flow. 

Using  the  knowledge  so  gained,  functions  representing  the 
velocity  components  were  constructed. 


II.  THE  VELOCITY  FIELD 

The  structure  of  the  velocity  field  was  graphically  shown 
by  visualization  techniques,  as  in  Figure  1.  It  is  this  pic¬ 
ture  which  clearly  demonstrates  the  vortical  nature  of  the 
flow  in  considerable  detail.  Starting  with  such  pictures  one 
can  develop  expressions  for  the  components  of  velocity,  ana¬ 
lytically  defining  the  flow. 

In  a  cross-section,  the  velocity  trajectories  are  a 
family  of  closed  curves.  Using  cylindrical  coordinates,  a 
conceivable  representation  of  these  trajectories  for  a  <p  — 
constant  plane  is: 
m=m^  n=n-^ 

k  =  \  )  a^  rm  zn  =  constant  (l) 

m=0  n=0 

In  non-dimensionalized  coordinates  a  particular  form  of  the 
series  is: 


A-.  Ap 

k  =  -(Z-R)  Z  “(1-R) 


A^  _  A A-,  > 

3  +  <Zc-Sc>  'Zo 
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which  is  suggested  by  the  equation  of  the  boundary. 


(Z-R)Z(l-R) 
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(3) 


where  the  non-dimensionalized  coordinates  and  velocities  are: 
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0  =  0 
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The  maximum  radius  of  the  cone  is  Tq  and  the  angular 
velocity  of  the  cone  is  «q  .  The  tangent  of  the  angle  be¬ 
tween  the  surface  of  the  cone  and  the  horizontal  surface  of 
the  cylindrical  container  is  h  .  In  the  future  this  angle 
will  be  referred  to  as  the  cone  angle  and  'will  be  denoted  by 
the  letter  a  .  Also  (R  ,  Z  )  are  the  coordinates  of  the 
vortex  center. 

For  these  trajectories  the  slope  of  any  tangent  line  is 
the  ratio  of  the  velocity  components  at  the  point  of  tangency. 
Also  the  velocity  components  should  be  such  as  to  satisfy  the 
continuity  equation.  Hence,  the  velocity  components  and 

u^  defined  in  terms  of  k  are: 


_  _E_  dk  _  5  3k 

R  c)Z  Z  R  SR 


where  E  is  a  function  of  k  ,  and  k  is  a  function  of 
and  Z  ,  given  by  equation  (2).  Thus, 


(b) 

R 
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E 

R 


(6) 


uR  =  -  ~  (i-r)*3z*2  1  (Z-R)^1  1[(X1+A2)Z  -  A2R] 


and 


Ur 


-  —  (l-R)^3  V2(Z-R)Xi  1[a3Z  -  (A1+  A3)R  +  Ax]  (7) 


Furthermore,  observation  shows  that  at  the  vortex  center,  u^ 
is  not  zero,  but  u^  and  u^  are  zero.  Hence,  from  the  ex¬ 
pressions  for  u^  and  u^  , 


(Ai  +  A2)Zc  -  A2Rc  =  0 


and 


(8) 


A3Zc  -  (An  +  AjR„  +  An  =  0 


3  c 


Consequently,  the  following  conditions  are  placed  on  A-^  , 


a2  >  a3 


K  ~  z„ 

c  c 


(9) 


Zc  -  Rc 

Re"1 


(10) 


Furthermore,  because  u^  and  u^  are  zero  on  the  boundary 
the  A ' s  must  be  chosen  so  that. 


Ax  -  1  >  0 


a2  -  1  >  0  , 


a3  -  1  >  0 


(11) 


It  follows  that  uR  will  be  zero  at  every  point  in  the  flow 
domain  where: 


Z  =  —  R 

R 

c 


(12) 


and  Ug  will  be  zero  at  every  point  where: 


Z  -  1  Z  -  R 
z  =  -2 - r - 2 - 2. 


Rc  -  X 


V1 


(13) 


For  values  of  Z  >  R  9  uR  is  positive  and  for  values  of 

c 

zc 

Z  <  ■£-  R  ,  Up  is  negative.  Likewise,  for  values  of 
nc 


Z  -  1  Z  -  R 

z  >  iihn  R  -  ~imr 

c  c 


]  ,  uz  is  positive  and  for  values  of 


Z  -  1  Z  -  R 

z  <  Hhn:  R  -  ttt 


]  ,  u7  is  negative.  These  results 


c  c 

agree  with  the  observations  using  flow  visualization  techniques. 


The  third  velocity  component  u^  could  conceivably  be 
represented  as  follows: 


u 


0 


=  Z(l-R) 


-f 


1-  z 


m=m^  n=n]_ 

Z  b  Djm 

m=0  n=0 


R*1  Zn 


(14) 


A  particular  simple  form  of  this  series  was  chosen  to  repre¬ 


sent  u. 

0 


as  follows: 


(15) 


u.  =  +  AZ(Z-R)  (l-R)  (r  +  bz) 

V  1  -z 

where  A  and  B  are  constants  to  be  determined.  Now  uR  , 
u0  ,  uz  have  been  determined  in  such  a  manner  that  the  con¬ 
tinuity  equation  and  the  boundary  conditions  are  satisfied. 

The  constants  ^  ^  B  and  the  parameter  E  can 

now  be  determined  by  experiment. 

III.  METHODS  FOR  THE  DETERMINATION 
OF  THE  Vs,  E,  A,  AND  B 

A.,  The  Determination  of  the  Vs 
The  Vs  can  be  determined  once  the  position  of  the  vor¬ 
tex  center  is  known.  The  following  is  a  description  of  the 
experimental  method  used  to  obtain  this  information. 

Small  water  particles  were  injected  into  the  flow  domain 
with  a  hypodermic  needle.  Using  two  micrometer  screws,  the  co¬ 
ordinates  of  the  particle  with  reference  to  the  top  and  side  of 
the  cylindrical  container  could  be  determined.  Water  particles 
moving  in  a  circular  path  about  the  axis  of  rotation  of  the 
cone  are  at  the  vortex  center.  For  castor  oil  U.S.P.  the  re¬ 
sults  are  presented  in  Table  I.  The  distance  of  the  vortex 
center  from  the  vertical  wall  is  £  ,  and  the  distance  of  the 
vortex  center  from  the  top  of  the  cylindrical  container  is  T|  . 
The  experimental  arrangement  for  the  cylindrical  container  and 
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one  of  the  cones  is  shown  in  Figure  2 


From  the  table  it  can  be  seen  that  as  the  cone  angle  de¬ 
creases,  the  vortex  center  moves  toward  the  cylindrical  wall. 
It  is  concluded  that  as  a  ,  the  cone  angle,  approaches  zero 
degrees,  R  approaches  unity.  As  a  increases,  the  vortex 

V 

center  moves  toward  the  axis  of  rotation  of  the  cone.  Some 
observations  made  with  a  76°  angle  cone  indicate  that  such 
is  the  case  also  for  cones  whose  angles  are  larger  than  45°. 

B.  The  Determination  of  E 
The  quantity  E  varies  with  the  surface  parameter  k 
and  the  angular  velocity  of  the  cone.  Its  determination  was 
made  from  considerations  of  what  may  be  called  types  of  cir¬ 
culation  of  a  particle  of  the  liquid.  Specifically,  the  cir¬ 
culations  considered  are  of  two  kinds;  one  around  a  cross- 
sectional  trajectory,  the  other  around  the  axis  of  rotation 
of  the  cone.  This  follows  from  the  fact  that  a  given  particle 
of  liquid  always  moves  on  a  particular  toroidal  surface.  In 
so  doing  the  particle  moves  not  only  around  the  fixed  circular 
axis  which  coincides  with  the  vortex  center  in  the  liquid  but 
also  around  the  axis  of  rotation  of  the  cone. 

The  time  required  for  a  particle  to  make  one  revolution 
about  the  fixed  circular  axis  is: 

t  -  £  SS.  (15) 
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where  ds  is  a  differential  length  on  a  particular  trajectory- 
defined  by  k  ,  and  v  is  the  magnitude  of  the  vector  sum  of 
the  velocity  components  v^  and  vz  .  Also,  the  above  time 
can  be  related  to  the  average  angular  velocity,  <oa^  ,  of  the 
particle  as  it  rotates  once  about  the  fixed  circular  axis. 

The  symbol  refers  to  a  particular  surface  defined  by  the 

parameter  k  .  The  relation  is: 

e  ’  /  *  (17) 


Non-dimensionalizing  the  integrand  and  solving  for  E  one 
finds  the  relation: 
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It  was  found  experimentally  that  the  average  angular 
velocity  of  the  particle  about  the  cone  axis.  Si  ^  ,  is  re¬ 
lated  to  the  average  angular  velocity  coa1  of  the  particle 
about  the  fixed  circular  axis  as  follows: 


co. 


ai 


=  R. 


si 


ai 


ai 


(19) 


where  Ra^  is  a  parameter  depending  only  upon  k  .  The  max¬ 
imum  angular  velocity  of  the  cone  for  the  velocity  range 


8- 


studied  is  0  .  Also  it  was  experimentally  found  that  n  . 

s  ai 

is  related  to  the  angular  velocity  of  the  cone  as  follows: 


n 


ai 


G  a 


0 


(20) 


where  G  is  a  parameter  depending  only  upon  k  .  Using  equa¬ 
tions  (19)  and  (20)  with  equation  (l8)  the  parameter  E  is 
given  as  follows: 


E 


Rnl(k)S(k)  Aq 
2t  Q 

s 


/ 


+  h2  dZ2] 
+  h2  uz2] 


1/2 

T72 


(21) 


The  two  parameters  Ra^  and  G  can  be  determined  experimen¬ 
tally  by  measuring  the  time  required  for  the  particle  to  travel 
once  about  a  trajectory  defined  by  k  ,  and  by  counting  the 
number  of  times  the  particle  moves  about  the  cone  axis  while 
moving  once  about  the  plane  trajectory.  The  integral  can  be 
evaluated  by  using  the  derived  velocity  functions. 

C.  The  Determination  of  A  and  B 
The  constants  A  and  B  can  be  calculated  using  data 
from  two  experimental  measurements.  The  first  measurement  is 
of  the  angular  velocity  of  water  particles  placed  at  the  vortex 
center  and  the  second  measurement  is  the  torque  on  the  cone 
surface. 


Using  a  stop  watch  over  long  periods  of  observation,  the 
angular  velocity  of  water  particles  placed  at  the  vortex  center 
was  determined.  The  results  are  shown  in  Table  I.  The  ratio 
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of  the  angular  velocity  of  the  particle  at  the  vortex  center 
to  the  angular  velocity  of  the  cone  is  .  The  values  in 
Table  I  are  the  averages  of  all  the  readings  taken  and  A ft ' 

V/ 

is  the  maximum  percent  of  variation  between  the  average  and 
an  individual  reading.  The  ratio  is  a  constant  over  a 

wide  range  of  angular  velocity.  The  range  for  which  the  re¬ 
sults  are  listed  in  Table  I  was  between  10  rpm  and  80  rpm. 

The  difference  in  individual  readings  is  probably  caused  by 
the  fact  that  the  water  particle  cannot  be  placed  at  the  exact 
center  of  the  vortex.  It  was  found  that  by  decreasing  the 
size  of  the  water  particle  as  the  angular  velocity  of  the 
cone  decreases,  the  agreement  between  individual  readings  was 
improved. 

As  the  cone  angle  approaches  zero  degrees,  the  angular 

velocity  at  the  vortex  center  approaches  one-half  the  value 

of  the  angular  velocity  of  the  cone.  As  the  angle  increases, 

the  ratio  o'  decreases.  For  a  =  45°,  the  ratio  is  0.306. 
c 

The  torque  on  the  rotating  shaft  was  measured  by  a  trans¬ 
ducer  constructed  and  developed  for  the  R.P.I.  Rotational  Fluid 
Flow  Generator.  The  description  of  this  transducer  and  the  co¬ 
related  electrical  circuit  is  completely  described  in  the  previ- 

2 

ously  mentioned  technical  report  .  The  torque  was  measured  for 
several  cone  angles  over  the  range  of  angular  velocity  from 
zero  to  100  rpm.  The  liquids  investigated  were  castor  oil 
U.S.P.,  S.A.E.  10W,  S.A.E.  20-20W,  and  S.A.E.  30  Lubrite  Motor 
Oil  manufactured  by  the  Mobile  Oil  Company.  Using  the  simplest 
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means  of  thermal  shielding  and  insulation,  the  experiments  were 
carried  out  at  a  temperature  of  22.2°  Centigrade.  The  results 
of  these  experiments  show  that  the  torque  is  a  linear  function 
of  the  angular  velocity  of  the  cone.  Hence,  the  constitutive 
equation  for  a  Newtonian,  incompressible  liquid  was  considered 
to  be  the  constitutive  equation  characterizing  the  liquids  for 
which  the  developed  velocity  field  is  applicable.  The  equation 
is  written: 

T  =  -  pi  +  2 p.X  (22) 

where  T  and  X  are  the  stress  tensor  and  the  rate  of  deforma¬ 
tion  tensor  respectively,  I  is  a  unit  tensor,  and  p  is  de¬ 
fined  as  the  arithmetical  mean  of  the  normal  stresses.  The  co¬ 
efficient  of  viscosity  is  p,  . 

In  cylindrical  coordinates,  the  shear  stress  component 
which  produces  a  torque  on  the  surface  of  the  cone  is  as  follows 

[Sv,  v,  Sv.  -i 

(_!  _  _l)Sin  p  -  — ^sinpj  (23) 

Sr  r  Sz  -I 

The  total  torque  on  the  surface  of  the  cone  is  then: 
a 

M  =  2tt  P  r^  [  Tq,  ]  •  sec  a  dr  (24) 

J  q  **  z=r  tan  a 

Putting  the  velocity  functions  into  equation  (24)  the  following 
formula  was  derived. 
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a 


iii=s£  -  2  ( i-R)  +  iog(i-R)-  is!  +  is!; 
I  2  5  6  ) 

(25) 

B  =  A(l  +  B)  (26) 

For  the  flow  examined  experimentally,  the  length  of  the 
free  surface  was  kept  at  10$  of  the  total  length  of  the 
largest  radius  of  the  cone.  Hence  in  equation  (25)  the  limit 
a  is  taken  to  be  0.9.  Then  equation  (25)  becomes: 

M/hn  •  cos  a  sin  a 

- v - - -  =  _  (.998  +  .0306)  (27) 

2TTiir0D 

Using  the  results  of  the  torque  experiments,  the  expression 
on  the  left  hand  side  of  equation  (27)  was  evaluated.  The  re¬ 
sults  are  shown  in  Table  II.  Using  the  average  value  for  each 
column,  6  can  be  calculated  for  each  angle  using  equation 
(27).  The  results  are  given  in  Table  III. 


M/£2q  •  cos  a  sin  a 
2^rn^ 


where : 


Using  the  values  of  6  in  Table  III,  the  values  for  the 
coordinates  of  the  vortex  center,  and  the  values  of  the  angular 
velocity  of  particles  placed  at  the  vortex  center,  given  in 
Table  I,  the  values  of  A  and  B  can  be  determined  as  follows 


=  Rc 

=  z_ 


Z„ ( 1-R„ )  AZ 

=  a  '  =  — - —  +  — £  (z  -R  )  ( l-R  )  (R  +BZ  )  (28) 

C  Rc(l-Zc)  Rc  «  cM  ° 
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(29) 


and. 


A  = 

1  +  B 

Solving  equations  (28)  and  (29)  simultaneously  A  and  B  are 
given  by  the  following  equations: 


A  -  - 


ft'  R„  ( 1-Z  ) 
c  c  c 


Z  fl-R  )  -  5  Z„  R  (Z  -Rj(l-R  )(1-Z  ) 
c  c  c  cxc  cv  CN  c 


z„  r  ( i-r )  ( l-z  )  ( z  -R  )‘ 

c  c  c  c  c  c 


(30) 


K  -  °  z  R02(zc-rc)(i-rc)(i-zc) 

°0  R0<1-z'0)  -  -  6  R0(zo-RcHi-Rc)  ( X-Z0 ) 

(31) 


IV.  ILLUSTRATION  OF  THE  METHOD  FOR  DETERMINATION 
OF  THE  VARIOUS  PARAMETERS 

Using  the  experimental  data  of  Table  I,  and  the  equations 
(9)  and  (10)  for  a  =  35°  ,  the  A's  are  as  follows: 

\  =  A2  =  4  *3  «  2  (32) 

The  equation  of  the  velocity  trajectories  is: 

k  =  -  (Z-R)V*(1-R)2  +  2.62  x  10-5  (33) 

In  Figure  3>  several  curves  have  been  plotted  using  fractional 
parts  of  2.62  x  10“^  for  k  .  A  comparison  can  be  made 
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between  the  calculated  trajectories  shown  In  Figure  3*  and 
the  actual  trajectories  shown  In  Figure  1. 

Using  the  experimental  data  given  In  Figures  4,  5,  6,  and 
the  curves  drawn  In  Figure  3*  the  parameter  E  was  determined. 
For  curve  "a",  E  =  13.6fiQ  ,  for  curve  "c",  E  =  8.65&Q  , 
and  for  curve  "d",  E  =  1.25^q  • 

Using  the  experimental  data  of  Tables  I,  II,  and  III,  for 
a  equal  35°  it  was  found  that: 

6  =  5.00 
K  "  *35 

R  =  .8  and  Z  =  .4  . 

c  c 

Using  equations  (30)  and  (31)  it  was  found  that: 

A  =  **19.6  and  B  =  *  1.25  . 

Hence,  the  velocity  components  are: 

uR  =  -iS  (1-R)2Z3(Z-R)3(2Z-R) 
n  R 

U„  =  -  £S  (l*.R)Z^(Z*R)3(Z-3R+2) 

Z  R 

v  =  •  19  .6Z(Z-R)  (l— R)  (Rwl.25Z) 

9  1  -z 


(34) 

(35) 

(36) 
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A  formula  for  finding  the  coefficient  of  viscosity  of 
liquids  using  a  35°  angle  cone  was  developed  using  the  re¬ 
sults  of  Table  II.  It  Is  as  follows: 

M,  =  .0256  eentipoise  (37) 

n0 

where  M  Is  in  dyne-cm  and  Is  In  secrl  .  The  coef¬ 

ficients  of  viscosity  calculated  from  the  formula  are  given 
In  Table  IV.  The  value  of  the  coefficient  of  viscosity  for 

castor  oil  U.S.P.  compares  well  with  that  given  In  the 

•3 

Smithsonian  Tables^  and  the  viscosities  of  the  S.A.E,  oils 

2i 

agree  with  the  values  given  by  Shaw  . 

Also  it  was  found  that  for  a  35°  angle  cone  the  circu¬ 
lation  time  about  the  fixed  axis  through  the  vortex  center 
decreased  as  the  vortex  center  was  approached.  For  £1q  equal 
to  40  rpm  the  circulation  time  for  curve  "a"  was  1.57 
minutes,  and  for  curve  "d"  the  value  was  .39  minutes.  An 
estimate  was  made  of  the  time  required  for  a  particle  to 
travel  about  a  curve  from  R  =  .3  to  R  =  .9  for  a  cone 

whose  angle  was  4°  .  For  equal  to  40  rpm  the  circu- 

26 

latlon  time  was  of  the  order  of  10  minutes.  Hence  as  the 
angle  of  the  cone  decreases  to  zero  It  is  concluded  that  the 
circulation  approaches  zero. 
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V.  EQUATIONS  OF  MOTION 

Applying  the  constitutive  equation  (22)  to  Newton's  equa¬ 
tions  of  motion  for  axisymmetric,  steady-state  flow  the  equa¬ 
tions  of  motion  which  have  the  linear  Stokesian  hypothesis  for 
a  constitutive  equation  are  the  Navier-Stokes  equations.  In 
cylindrical  coordinates,  the  non-dimensionalized  equations  for 
the  specific  problem  under  consideration,  neglecting  body 
forces  are: 


u 


R 


dUp 

dR 


+  ur 


^R 

dZ 


1  f  dP 
K  I-  dR 


d2u 


+ 


R 


dZc 


*■» 

dR2 


1  ^R  _  ur1| 
R  dR  I?" 


(38) 


u 


R 


^0 

dR 


+ 


URU0 

R 


i  du . 
+  & 
R  dR 


where  P  is  related  to  the  arithmetical  mean  p  of  the  normal 
stresses  by  the  relation 


f  S  ^0P 

p(r,  z)  =  — 5- 

n 


and  K  is  defined  as. 


-1 6- 


(41) 


The  Reynolds 1  number  Re  Is  proportional  to  K  and 
given  as  follows; 

Re  =  (43) 

h 

The  Reynolds 1  number  refers  only  to  a  cone  of  given  angle  and 
given  gap  width  between  the  cone  and  cylindrical  wall  at  the 
free  surface.  For  a  35°  angled  cone,  the  corresponding 
range  of  the  Reynolds 1  number  using  castor  oil  Is 

.48  <  Re  ^  48  (44) 

The  theoretical  velocity  field  closely  approximates  the 
observed  velocity  trajectories  in  a  <f>  =  constant  plane,  and 
satisfies  the  continuity  equation  and  the  boundary  conditions. 
Also,  the  derived  velocity  components  Indicate  that  the  parti¬ 
cles  of  the  fluid  spiral  around  the  central  curve  of  a  vortex 
as  they  move  about  the  axis  of  the  cone.  Also,  the  theoretical 
velocity  components  agree  with  the  experimentally  determined 
values  at  the  vortex  center.  Using  the  velocity  component  in 
the  0  direction,  the  total  torque  calculated  on  the  cone  is 
approximately  the  same  as  that  calculated  on  the  cylindrical 
container.  It  is  on  the  basis  of  the  agreement  of  these  re¬ 
sults  with  the  actual  flow  produced  in  the  laboratory  that  the 


theoretically  determined  velocity  field  is  considered  to  be 
an  approximate  solution  of  the  differential  equations  ( 38 ) > 
(39)  and  (40). 

Furthermore,  integrating  equation  (40)  with  respect  to  Z 
and  equation  (38)  with  respect  to  R  ,  two  representations  of 
the  non*dimensional  pressure  field  are  obtained. 


is 

dR 


-2-  +  uz  —  TV  dR  +  F,(Z) 

R  Z  dZ  ,J  1 


(45) 


If  the  velocity  components  u^  »  uR  >  uz  are  exa°t  solu¬ 
tions  of  the  Navi er- Stokes  equations  then  should,  of 

course,  equal  P2  • 

The  discrepancy  that  may  exist  between  P1  and  P2  in 
the  present  investigation  has  not  been  examined.  It  is  an  ob¬ 
ject  of  study  for  future  research. 
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Finally  it  is  suggested  that  an  infinite  series  expansion 


of  k  and  u,  of  the  form; 

0 


k  = 


00  00 


n^t>  n^b 


C  Rm  Zn 
mn 


(47) 


and 


00  00 


U  . 


Z(l»R) 


U *£W  )  ) 

i  -  z  _f-u  M 


D  Rm  Zn 
mn 


(48) 


m=0  n=0 


may  provide  an  exact  solution  to  the  equations  of  motion. 


VI.  DISCUSSION  AND  CONCLUSION 
Using  observations  of  the  actual  flow  patterns  and  measure¬ 
ments  of  torque  and  angular  velocity,  it  was  possible  in  the 
present  study  to  determine  satisfactorily  the  functions  describ¬ 
ing  the  velocity  field  induced  in  a  liquid  by  a  rotating  cone. 
The  stress  field  and  the  constitutive  equations  were  studied  in 
terms  of  the  velocity  functions.  The  relation  between  the  angu¬ 
lar  velocity  of  the  cone  and  the  total  torque  was  shown  to  be 
linear.  Therefore,  it  was  concluded  that  the  liquids  used  in 
the  investigation  were  Newtonian  and  the  differential  equations 
of  motion  governing  the  flow  are  the  classical  Navier-»Stokes 
equations.  Because  of  the  many  agreements  of  the  constructed 
velocity  functions  with  the  actual  flow  for  cones  of  any  angle. 
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they  are  considered  to  be  reasonable  approximations  to  those 
which  are  a  solution  to  the  Navier-Stokes  equations.  As  the 
angle  of  the  cone  approaches  zero  as  a  limit,  the  circulation 
about  the  vortex  center  was  shown  to  approach  zero.  The  torque 
and  angular  velocity  relation  approaches  that  which  has  been 
developed  by  rheologists  for  the  determination  of  the  coef¬ 
ficient  of  viscosity  for  Newtonian  liquids,  when  rotating  the 
liquids  by  cones  whose  angles  are  less  than  4°  .  However,  by 
using  the  velocity  functions  developed  in  the  present  research, 
the  measurement  of  viscosity  is  now  not  limited  to  the  use  of 
cones  with  small  angles  and  small  amounts  of  liquid  but  can  be 
carried  out  with  cones  of  any  angle  rotating  in  any  amount  of 
liquid. 

It  Is  finally  concluded  that  the  investigation  provides 
effective  new  knowledge  of  the  flow  generated  by  a  rotating 
cone  .immersed  in  a  liquid  and  furthermore,  provides  the  rational 
basis  for  the  generalization  of  formulae  now  used  in  visco- 
simetry  for  Newtonian  liquids. 
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Figure  1.  Structure  of  the  Flow 
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Figure  5-  PLOT  OF  R  .  AS  A  FUNCTION  OF 
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Measurement  of  Coordinates  of  the  Vortex  Center  and  of 
Angular  Velocity  of  Particles  at  the  Vortex  Center. 


o 

in 

on 

ii 

vd 

o 

• 

o 

00 

rH 

rH 

1 

on 

o 

+ 

00 

on 

rH 

l 

ON 

o 

• 

ON 

rH 

1 

■ 

El 

EB 

IRfl 

■ 

CM 

o 

• 

CM 

rH 

CM 

1 

■ 

CM 

OH 

• 

CM 

■ 

CM 

a 

lO 

c 

<D 

U 

O 

<M 


U 

o 


CQ 

0) 

G 

o 

o 


aO  0) 


<4H 

O 

« 

0> 

2 

rH 

S> 


tf) 

5 


Table  III 


T  C  Material 


DISTRIBUTION  LIST 


Copies 

65 

10 


6 

6 

3 

1 

1 

1 

2 

1 

2 

1 

1 

1 


Organization 

Commanding  Officer  and  Director,  David  Taylor  Model 
Basin,  Washington,  D.  C.  20007  Attn:  Code  513 

Chief,  Bureau  of  Ships,  Department  of  the  Navy, 
Washington  25,  D.  C. 

3  -  Technical  Library 
1  -  Preliminary  Design  (Code  420) 

1  -  Hull  Design  (Code  440) 

1  -  Research  and  Development  Planning  Division 
(Code  330) 

1  -  Hydromechanics,  Logistics  &  Special  Craft 
(Code  341B) 

1  -  Ship  Silencing  (Code  345) 

1  -  Laboratory  Management  Division  (Code  320) 

Chief,  Bureau  of  Yards  and  Docks,  Department  of  the 
Navy,  Washington,  D.  C. 

Chief,  Bureau  of  Naval  Weapons,  Department  of  the 
Navy,  Washington  25,  D.  C. 

3  -  RU 
3  -  RAAD-3 

Chief  of  Naval  Research,  Department  of  the  Navy, 
Washington  25,  D,  C.  Attn:  Fluid  Dynamics  (Code  438) 

Director,  U.  S.  Naval  Research  Laboratory, 

Washington  25,  D.  C.  Attn:  Code  2021 

Commanding  Officer,  Office  of  Naval  Research  Branch 
Office,  230  N.  Michigan  Avenue,  Chicago  1,  Illinois 

Commander,  U.  S.  Naval  Ordnance  Laboratory,  White  Oak, 
Silver  Spring,  Maryland 

Commander,  U.  S.  Naval  Ordnance  Test  Station,  3202  East 
Foothill  Boulevard,  Pasadena,  California 
1  -  Dr.  A.  G.  Fabula 

Commanding  Officer  and  Director,  U.  S.  Navy  Underwater 
Sound  Laboratory,  Fort  Trumbull,  New  London,  Conn. 

Director,  Ordnance  Research  Laboratory,  Pennsylvania 
State  University,  University  Park,  Pennsylvania 

Superintendent,  U.  S.  Naval  Postgraduate  School, 

Monterey,  California  Attn:  Librarian 

Assistant  Secretary  of  the  Army,  Research  and  Develop¬ 
ment,  Washington  25,  D,  C. 

Director,  U.  S.  Waterways  Experiment  Station,  Corps  of 
Engineers,  P.O.Box  631,  Vicksburg,  Miss., Attn:  Tech.Libn. 


Organization 


Office  of  the  Chief  of  Engineers,  Civil  Works, 

Department  of  the  Army,  Oravelly  Point, Washington  25,D.C 

Director  of  Research,  National  Aeronautics  &  Space 
Admin.,  1525  H  St.,  N.W.,  Washington  25,  D.  C. 

Director,  Langley  Aeronautical  Laboratory,  National 
Aeronautics  &  Space  Admin.,  Langley  Field,  Virginia 

Director,  Hydraulic  Laboratory,  Bureau  of  Reclamation, 
Denver  Federal  Center,  Denver,  Colorado 

Defense  Documentation  Center-  for  Scientific  &  Technical 
Information,  Arlington  Hall  Station,  Arlington  12,  Va. 

Director,  National  Bureau  of  Standards,  National 
Hydraulic  Laboratory,  Washington  25,  D.  C. 

Newport  News  Shipbuilding  &  Dry  Dock  Co.,  Newport  News, 
Virginia 

1  -  Assistant  Naval  Architect 
1  -  Director,  Hydraulic  Laboratory 

California  Institute  of  Technology,  Division  of 
Engineering,  Pasadena  4,  California 

Director,  Davidson  Laboratory,  Stevens  Institute  of 
Technology,  711  Hudson  St.,  Hoboken,  New  Jersey 

Director,  Hydraulic  Laboratory,  Worcester  Polytechnic 
Institute,  Worcester  2,  Massachusetts 

Director,  Hydrodynamics  Laboratory,  Massachusetts 
Institute  of  Technology,  Cambridge  39»  Massachusetts 

Head,  Dept,  of  Naval  Architecture  &  Marine  Engineering, 
Massachusetts  Institute  of  Technology,  Cambridge  39> 
Massachusetts 

Director,  Woods  Hole  Oceanographic  Institute, 

Woods  Hole,  Massachusetts 

Western  Co.  of  North  America,  1171  Empire  Central, 

Dallas  7>  Texas  Attn:  Dr.  H.  R.  Crawford 

Director,  Engineering  Societies  Library,  29  West  39th 
Street,  New  York  18,  New  York 

Library,  California  Institute  of  Technology, 

Pasadena,  California 

Librarian,  Massachusetts  Institute  of  Technology, 
Cambridge  39,  Massachusetts 


Copies 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

2 


2 

1 

1 

1 

1 

1 


( 


Organization 

Librarian,  School  of  Engineering,  University  of 
Texas,  Austin,  Texas 

Serials  Division,  University  of  Minnesota  Library, 
Minneapolis,  Minn. 

Prof.  A.  B.  Metzner,  Dept,  of  Chemical  Engineering, 
University  of  Delaware,  Newark,  Delaware 


Prof.  P.  H.  Clauser,  Dept,  of  Mechanics, 
The  Johns  Hopkins  Univ.,  Baltimore  18,  Md. 


Prof.  S.  Corrsin,  Dept. 
The  Johns  Hopkins  Univ,,, 


of  Mechanics, 
Baltimore  18, 


Md. 


Dr.  P.  Kaplan,  Oceanics,  Inc.,  Technical  Industrial  Pk. 
Plainview,  Long  Island,  New  York 


Dr.  John  We3ke,  Hydro.  Div. 

Univ.  of  Maryland,  College  Park,  Maryland 


Dr.  George  A.  Wislicenus,  Navy  Ord.  Water  Tunnel, 
Penn  State  Univ.,  State  College,  Pa. 

Dr.  A.  H.  Shapiro,  M.I.T. 

Cambridge,  Mass . 

Dir.  Iowa  Inst,  of  Hydraulic  Res. 

State  Univ.  of  Iowa,  Iowa  City,  Iowa 
Dr.  L.  Landweber 


Dr.  Karl  E.  Schoenherr 

7053  Western  Ave.,  Washington,  D.  C. 

Dept,  of  Eng.,  Inst,  of  Eng.  Res., 

Univ.  of  California,  Berkeley,  California 
1  -  Dr.  J.  Wehausen 

Hydronautics,  Inc.,  Pindell  School  Road, 

Laurel,  Howard  County,  Maryland 

Dr.  J.  L.  Erlcksen 

The  Johns  Hopkins  Univ.,  Baltimore  18,  Md. 

Dr.  W.  Philippoff,  Esso  Res.  &  Eng.  Co., 

Linden,  N.  J. 

Prof.  E.  T.  Chu,  Dept,  of  Eng. 

Yale  Univ.,  New  Haven,  Conn. 

Dr.  H.  Markovitz,  Mellon  Inst. 

Pittsburgh,  Pa. 

Prof.  John  D.  Perry,  Univ.  of  Wisconsin,  Madison,  WIs. 


